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We measure the upper state lifetime and two ratios of vibrational branching fractions f v i v on 
the B 3 ni(u') - X 1 S+(u) transition of TIF. We find the B state lifetime to be 99(9) ns. We also 
determine that the off-diagonal vibrational decays are highly suppressed: /01//00 < 2 x 10 -4 and 
/02//00 = 1.10(6)% , in excellent agreement with their predicted values of /01//00 < 8 x 10~ 4 and 
/02//00 = 1.0(2)% based on Franck-Condon factors calculated using Morse and RKR potentials. 
The implications of these results for the possible laser cooling of TIF and fundamental symmetries 
experiments are discussed. 



The laser cooling of molecules presents a daunting chal- 
lenge with potentially rich rewards. Transverse cooling 
and collimation of a cold molecular beam can be accom- 
plished by scattering a few hundred photons from each 
molecule, while over 10,000 photons must be scattered 
to bring a typical molecule in the beam to rest. It is 
difficult to find molecular cycling transitions that will 
allow so many absorption and decay cycles without sig- 
nificant loss to the myriad of rotational and vibrational 
states present in most molecular systems. In addition, 
successful laser cooling requires that the excited state be 
relatively short lived, such that it can complete many 
absorption/emission cycles before leaving the laser inter- 
action region. Despite these challenges, transverse cool- 
ing has recently been achieved in strontium monoffiioride 

In the present paper we explore the possibility of 
laser cooling thallium monofluoride (TIF). High precision 
searches for the Schiff moment and the proton electric 
dipole moment (EDM) have been carried out in beams 
of TIF d, 3] ■ These experiments are tests of both par- 
ity (P) and time-reversal (T) symmetries. TIF exhibits 
a large enhancement in its sensitivity to such violations 
due to the large internal electric field of the molecule 
and thallium's large atomic number Z = 81 []|. It also 
displays a remarkable insensitivity to systematic effects 
associated with external magnetic fields. The TIF exper- 
iments were largely limited by the relatively broad line 
widths associated with the molecules' rapid transit time 
through the apparatus, and the modest thermal popu- 
lations of the state of interest. Using cryogenic beams 
and laser cooling, it may be possible to overcome these 
limitations and TIF might again emerge as an interest- 
ing candidate for measuring symmetry violations in the 
nucleus. 

The transition = 0, J p = 1~) - B 3 Ui(v' = 

0,J' P = 1 + ) (where J p denotes the rotational angular 
momentum and parity) of TIF is an interesting candi- 
date for a cycling transition. As we argue here, it ap- 
pears that this transition should be highly closed to other 
(unwanted) electronic, vibrational, and rotational decay 
paths. We begin with a discussion of electronic decay 



paths. The only other electronic decay transition from 
the B 3 Hi state is to the A 3 n + state (and its as-yet un- 
observed 3 n - partner). The branching fraction for this 
transition should be very small, according to the follow- 
ing logic. First, the transition frequency ujba for the B-A 
transition is much smaller than that for the B-X tran- 
sition (ljbx)- ^bx/ojba ~ 22. Moreover, the electric 
dipole (El) matrix element dsA should be several times 
smaller than dsx- This can be seen as follows. Nomi- 
nally, both the B-X and B-A El transitions are forbid- 
den: the former requires a change in the total spin S and 
the latter a change in the projection £ of S along the in- 
ternuclear axis. Both transitions acquire a non-zero ma- 
trix element due to off-diagonal mixing with other elec- 
tronic states, via the spin-orbit (SO) interaction. For the 
B-X transition, this occurs primarily due to mixing of the 
nearby 6sa6pir C 1 Hi state into the QsaQpn B 3 Iii state; 
the admixed QsaQpir C 1 Hi — (6scr) 2 X 1 I] + Rydberg tran- 
sition is very strong. However, due to selection rules for 
SO mixing (AS = ±1,0; AS = ±1,0; AA = ±1,0 , where 
A is the projection of orbital angular momentum along 
the internuclear axis; and Af2 = 0, where f2 = A + £) [6], it 
can be verified that the B-A transition requires SO mix- 
ing with another orbital configuration (or a second-order 
mixing if only the lowest-lying configurations are consid- 
ered) . The nearest relevant level is likely the as yet unob- 
served 6sa6pa 3 £ + state, which can mix into both the B 
state and the 3 n - state; however, this admixture leads 
only to a weak perpendicular-band 0] (6p7r — 6po~) tran- 
sition. Taking into account all possibilities, we crudely 
estimate dux I Aba ~ 3. Finally, since the decay rate 
r oc d 2 uj 3 , we estimate that Tbx/^ba ~ 10 5 . As such, 
unwanted B-A decays are unlikely to significantly limit 
the cooling process. 

Selection rules limit the rotational states of X 1 T, + ac- 



cessible to the decay from the B 3 Hi( 
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state. In the absence of hyperfine structure (HFS), this 
state can decay only to X 1 E + (J P = 1~). Off-diagonal 
mixing of different rotational levels due to HFS can in 
principle allow unwanted rotational transitions. In the 
X 1 ^ state the HFS is very small (~kHz)Q and such 
mixing is negligible. In the B 3 Xl\ state, the HFS in- 



2 



teraction can be estimated from a simple scaling law 
(HFS is roughly proportional to Z) combined with the 
observed HFS in the isoelectronic species A1F H . Based 
on this, HFS interaction with the Tl nucleus (Jti = i 
for both 205 T1 and 203 T1) can be substantial, with es- 
timated B 3 Ui(v' = 0, J' p = 1+) state HFS interaction 
strength of order At\ ~ 1 GHz, while HFS associated 
with the 19 F nucleus (If = |) should be much smaller 
(Ap » 10 MHz). Hence we dehne the intermediate an- 
gular momentum Fi = J + iri and the total angular mo- 
mentum F = F! +I F . The B 3 Ui(v' = : J' p = 1+) state 
should split into well-resolved levels with F{ = i and | , 
each of which will consist of a pair of closely-spaced lev- 
els with F' = F[ ± \. Cycling on the X 1 S+(J P = 
l-,Fi,F) - B 3 Ili(v' = 0,J' p = l + ,F{ = i,F') manifold 
of nearly-degenerate transitions should have minimal loss 
to other X 1 E + (J") levels, since these B state sublevels 
mix with levels of higher J' only due to the weak HFS 
interaction of 19 F. In this case, the branching fraction 
fj'j" to unwanted higher rotational states can be esti- 
mated as fjijn ~ (A F /Sf ) 2 < 1CT 6 , where Bf is the 
rotational constant in the B^Iii state. Thus, using this 
particular HFS transition should enable cooling with lit- 
tle loss. 

Because no selection rules limit the various possible 
vibrational transitions, branching to unwanted vibra- 
tional levels is typically the most challenging problem 
in the laser cooling of a molecule @. We have numer- 
ically calculated Franck-Condon factors (FCF) for the 
low lying vibrational transitions under various assump- 
tions and approximations, in an attempt to quantify the 
uncertainties in their values. In all cases, precise spec- 
troscopic data, suitable for determining both the £? 3 Hi 
and X 1 Tj + state potentials, is taken from Ref. [l(|. In 
one approach, we used the full set of tabulated Dunham 
coefficients for both states, using standard routines for 
determining RKR potentials and their associated FCFs 
[ill . However, we believe there are two significant am- 
biguities in this approach. First, the difference in effec- 
tive rotational constants between e and /-parity levels in 
the B 3 Xl\ state, due to ^-doubling, has a non-negligible 
effect on the FCFs at the precision of interest here; it 
is not clear whether this difference is physically mean- 
ingful in terms of its effect on the molecular potentials. 
Second, the data in Ref. [l(| includes the Dunham coef- 
ficient Y30 (second vibrationally anharmonic term) only 
for the B 3 ni state, while the Dunham expansion termi- 
nates at Y20 for the X 1 Y, + state; it is unclear whether 
this inherent asymmetry in treatment of the potentials 
can lead to errors in the FCFs, particularly for higher 
vibrational levels of the X state. Hence, we evaluated 
FCFs for RKR potentials with rotational Dunham co- 
efficiencts of both e and /-parity levels, and both with 
and without inclusion of the Y30 term for the B state. 
Finally, in order to quantify possible errors in the FCFs 
due to the stated uncertainties in the spectroscopic data, 
we evaluated the FCFs using Morse potentials and the 
associated analytic form of the vibrational wavefuctions 



[l2l |. Here we used the Dunham coefficients Y10 and Y20 
for the harmonic and anharmonic vibrational terms of 
the Morse potential; for the X state we used the Dun- 
ham coefficient Ybi as the effective rotational constant to 
determine the internuclear separation r e — 2.084438 A, 
while for the B state we used r' e = 2.0740(5) A as de- 
rived in Ref. [l(| from a combined potential fit to both 
e- and /-parity levels. Within this model we calculated 
FCFs with all input parameters varied within their stated 
ranges of uncertainty. 

We quote a total uncertainty range for our calculated 
FCFs that incorporates all values obtained from these 
various calculations. The resulting Franck-Condon fac- 
tors are shown in TablelH The FCF matrix is highly diag- 
onal, so it appears that the X(v = 0) to B(v' — 0) transi- 
tion could provide a good cycling transition, with leakage 
to other vibrational levels broadly similar to that of the 
demonstrated case of SrF 0,0 ■ However, the uncertain- 
ties in the FCFs are significant at the level needed to eval- 
uate exactly how many, and which, vibrational repump- 
ing lasers would be needed in order to scatter enough pho- 
tons for laser cooling and/or slowing of TIF. In addition, 
the B 3 Hi potential of TIF, while mostly arising from an 
ionic bond between Tl + and F~ , is substantially modi- 
fied by a curve crossing with a covalent-bonding potential 
[l3|]. Hence, it seems conceivable that the X-B transition 
dipole moment could change significantly with the inter- 
nuclear distance; in this case the FCFs alone would not 
be sufficient to determine the vibrational branching frac- 
tions. Given all the possible sources of error, we chose to 
measure two of the vibrational branching ratios to test 
the accuracy of our calculations. 

Our measurements were made using a molecular beam 
that originates from a stainless steel oven filled with TIF. 
The oven is heated to temperatures in the range 415-460 
C, vaporizing the TIF which escapes the oven through a 
hole in the chamber. A set of four hollow ceramic tubes 
(25 mm long, 2.4 mm id and 3.3 mm od) fill this opening. 
These tubes are kept warmer than the oven and serve to 
pre-collimate the beam. The beam is further collimated 
by a 6 mm high by 10 mm wide aperture located about 
30 cm from the oven. The interaction region where the 
TIF beam intersects our laser beam is located about 6.5 
cm downstream from this aperture. This entire beam 
assembly is contained in a cylindrical stainless steel vac- 
uum chamber which is maintained at a pressure of about 
10~ 6 Torr by diffusion pumps. 

In order to excite the X(v = 0) to B(v' — 0) transition 
we need to generate ultraviolet (uv) light at 271.7 nm. 
We use a tripled YAG laser (Quantel Model 770B) at 355 
nm to pump a dye laser (Quantel Model TDL60) with 
Coumarin 540A dye. The output of this laser is tunable 
near 543.4 nm. The beam is directed to a BBO crys- 
tal where it is frequency doubled to produce the desired 
uv wavelength. The residual green light emerging from 
the BBO crystal is directed to a wavemeter (New Focus 
Model 7711) where its wavelength is monitored. The uv 
laser beam is incident upon the molecular beam at a right 
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0.989(2) 


< 0.0008 


0.011(2) 


< 0.0003 


< 0.0002 


< 0.0003 


< 0.0002 
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< 0.0004 


0.949(6) 


0.016(8) 


0.031(7) 


< 0.006 


< 0.0011 


< 0.002 


2 


0.010(2) 


0.03(2) 


0.77(6) 


0.09(5) 


0.08(2) 


0.03(2) 


< 0.006 



TABLE I: Calculated Franck- Condon Factors for the X 1 E + — _B 3 ITi transition of TIF. The primed (unprimed) vibrational 
quantum numbers refer to the excited (ground) electronic level. As described in the text, the central values and the uncertainties 
have been chosen to accommodate the variations in the predictions from the different model potentials and the uncertainties 
in the spectroscopic parameters used as input to the models. 



Branching ratio 


Theory 


Experiment 


/01//00 


< 0.0008 


< 0.0002 


/02//00 


0.010(2) 


0.0110(6) 



TABLE II: A comparison of our measured and predicted ra- 
tios of branching fractions from the -B 3 IIi, v' = state to the 
various vibrational states v of the X E + level in TIF. The 
branching fraction is taken as proportional to the Franck- 
Condon factor times the cube of the decay energy for that 
branch. 



angle and their overlap volume defines the region of the 
molecular excitation. To minimize scattering of the laser 
light, the linearly polarized laser beam enters and exits 
the vacuum chamber through quartz Brewster windows, 
mounted about 38 cm from the molecular beam. 

Detectors are placed above and below this interac- 
tion region, on a line perpendicular to both the laser 
and molecular beams and parallel to the laser polariza- 
tion. Fluorescence from the B state passes out of the 
vacuum chamber through a window, is collimated by a 
lens, passes through an interference filter and is focused 
on an adjustable diameter aperture that serves as a spa- 
tial filter. The fluorescence which passes through the 
aperture is detected by a uv sensitive phototube. All 
optics are fused silica to allow transmission of the uv 
light. Following preamplification, the single photon sig- 
nals are counted by a photon counter (Stanford Research 
model SR400). The device counts the number of pho- 
tons collected by each of our two detectors and relays 
this information to a computer for storage. 

We use the lower detector to constantly monitor the 
strong fluorescence signal at 271.7 nm. This provides 
a useful relative calibration of the fluorescence intensity 
and allows us to remove fluctuations in laser intensity 
and frequency and molecular beam intensity from our 
measurements. To measure branching ratios, we alter- 
nate between three different interference filters in the 
upper optics assembly. One of these is chosen to trans- 
mit the strong v' — to v = fluorescence near 271.7 
nm. The other interference filters allow us to monitor 
respectively the v' = to v = 1 and v' = to v = 2 de- 
cays at 275.2 nm and 278.8 nm. In order of increasing 
wavelength these three, 2 inch diameter filters (Andover 
Corporation) have center wavelengths of 271.6 nm, 275.3 
nm, and 279.3 nm, and FWHM bandwidths of 10 nm, 



1.5 nm and 2.0 nm. We tune the laser to achieve the 
largest possible fluorescence signal. This occurs at 271.7 
nm in a region of the rotational spectrum where many 
of the low rotational number Q-branch transitions are si- 
multaneously excited by our broad-band laser. We refer 
to this region as the Q-branch bandhead. We repeatedly 
compare the numbers of photon detections with each of 
our three filters, with the excitation laser tuned both on 
and off resonance. The numbers are corrected for leakage 
of the strong transition at 271.7 nm through the other fil- 
ters. Taking into account the transmissions of the filters 
at the various wavelengths and the finite angular spread 
of the fluorescence passing through the filters, we infer 
the branching ratios reported in Table [TTJ The agree- 
ment with the predicted values is excellent. Figure Q] 
summarizes our calculated and experimental results for 
the branching fractions for the decay of the B 3 Hi(v' — 0) 
state. Nearly 99% of the excited state population should 
return to the v = level and about 1% will return to the 
v = 2 level. The decays to v — 1, 3, 4, 5 and 6 or greater 
should all be less than a few hundredths of a percent. 

The B state lifetime also plays a critical role in as- 
sessing possible laser cooling schemes. To measure the 
lifetime, we feed our photomultiplier output into a multi- 
channel scalar (Stanford Research SR430). The instru- 
ment sorts the emitted photon counts into sequential 5 
ns bins. We fit this data with three parameters describ- 
ing an exponential decay plus a constant background. 
Care is taken to avoid photon pile-up problems at the 
beginning of the decay. We measure the lifetime on the 
Q-branch bandhead as well as on the isolated rotational 
lines R23, R33, R43 and P61. Each measurement is done 
using two different phototube/amplifier assemblies. As a 
further check, additional lifetime measurements are made 
by simply averaging the decay fluorescence on a digital 
300 MHz oscilloscope for the Q-branch bandhead and the 
lines R13, R23 and R33. The statistical uncertainty in 
the lifetime associated with any individual fit (Figure [5]) 
is typically less than a few percent. We observe some vari- 
ations in the fit lifetimes with the detection and analysis 
system used and the selection of the time interval chosen 
for the fit. The majority of these variations appear to 
be associated with our inability to completely remove an 
electronics background associated with the firing of the 
laser Q switch. We have chosen to quote an uncertainty 
in our measured lifetime that encompasses all of the ob- 
served variations. Combining the results from all of our 
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FIG. 1: The possible decay paths of the v' = level of the 
B 3 Hi state. The values we have deduced for the branching 
fractions fo v to the various vibrational levels of the X E + 
state are shown. Our estimated upper bound on the branch- 
ing fraction (bf) to the j4 3 IT + state is also shown. 




200 300 400 

Time (ns) 



T" 

500 



600 



FIG. 2: (Color online.) A typical lifetime measurement on the 
Q-branch bandhead. The points are the number of photon 
counts recorded in successive 5 ns bins. The initial time in 
the fit and plot is chosen to be 110 ns after the laser pulse. 
The data has been accumulated over 32,767 laser pulses and 
corrected for an offset and an electronics background. The 
solid line is a simple exponential fit with a 101 ns decay time. 



measurements we conclude that the excited state lifetime 
is 99 (9) ns. To the best of our knowledge, this is the first 
lifetime measurement of the B state of TIF. 

Next we consider the implications of our results for 
laser cooling and slowing of a TIF beam, and what this 
would mean for future experiments using TIF in tests of 
fundamental symmetries. We discuss two possible sce- 
narios: one which simply uses an improved molecular 



beam, and another in which the TIF molecules could 
be slowed and trapped for much longer coherence time. 
First, consider the possibility of using lasers to trans- 
versely cool and collimate a beam of TIF. Our results 
suggest that with a cooling laser tuned to the Ql line 



X(v = 0, J=1,F!,F) to B(v' = 0, J 1 -- 



1,F>: 



,F') tran- 



sition at 271.7 nm and a single repump laser, tuned to 
the analogous HFS component of the Ql line but on the 
X(v = 2) to B(v' = 0) transition at 278.8 nm, one should 
be able to achieve several thousand cycles on the cooling 
transition. This will be sufficient for transverse cooling 
and collimation of a cryogenic TIF beam [l4j]. Such a 
beam would have significantly increased occupation of 
the low rotational and vibrational levels of the molecule, 
as well as significantly lower forward velocity, compared 
to previous work with TIF beams. Based on typical 
brightness, velocity distribution, and internal tempera- 
ture achieved with this type of molecular beam source 
14{ , plus a typical transverse velocity capture range (cor- 
responding to a transverse temperature ~ 5 mK) and 
final transverse temperature (< 300/iK) of laser cool- 
ing [l]], we estimate crudely that an EDM experiment 
conducted in such a collimated cryogenic beam might 
achieve approximately two orders of magnitude improve- 
ment over the present TIF limit. This would yield limits 
on the proton EDM and Schiff moments comparable to 
those that have been achieved in the Hg EDM experiment 
15]. 

Dramatic improvement in interaction time could be 
achieved if the beam can be stopped and trapped. A 
TIF cryogenic beam would likely have a mean veloc- 
ity of about v w 150 m/s 14]. Approximately N — 
mvX/h w 23,000 photon absorption/emission cycles are 
needed to stop these molecules. To achieve this large 
number of cycles will likely require additional repump 
lasers to keep the molecules from accumulating in the 
X(v= 1,3,4 and/or 5) levels. These levels can be re- 
pumped either through excitations to the B(v' = or 2) 
levels or the A 3 n 0+ (v' = 0) state (see Table UnJ) [ll|. 
Hence, with a few additional lasers it is likely that the 
beam could be stopped using the radiative force from 
photon scattering. Further measurements of the vibra- 
tional branching fractions for decays from B(v' — 0) to 
AT(u = l,3,4 and/or 5), at a level of sensitivity beyond 
what we were able to achieve, will be needed to evaluate 
exactly what lasers will be needed for such a task. 

With knowledge of the excited state lifetime, we can 
estimate the length of the apparatus required to carry 
out radiative force slowing. For a two level system, op- 
timized spontaneous cooling allows one photon cycle to 
be completed in every two lifetimes of the excited state. 
In our X(J = 1) level of the ground state, the exis- 
tence of dark Zeeman states will cause the photon cy- 
cling rate to decrease by about a factor of three Q. 
(The remixing of these dark states into the optical cy- 
cle can be accomplished by rapid switching of the slow- 
ing lasers' polarization |17| or by resonant microwave 
transfer via the X(J = 0) level |2|.) We estimate the 
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0.89(1) 


0.087(1) 


0.018(1) 


0.0024(3) 


0.0003(1) 


0.00004(1) 


1 


0.11(1) 


0.75(2) 


0.10(1) 


0.036(1) 


0.005(1) 


0.0010(2) 


2 


0.0003(1) 


0.15(1) 


0.71(2) 


0.075(6) 


0.051(3) 


0.007(1) 



TABLE III: The Franck-Condon Factors for the - ,4 3 n o+ transition of TIF calculated using Morse potentials derived 

from the spectroscopic parameters of Ref. The primed (unprimed) vibrational quantum numbers refer to the excited 

(ground) electronic level. The uncertainties in the last digits listed in parentheses reflect the uncertainties in the spectroscopic 
data used to determine the Morse potential and do not reflect any limitations associated with the Morse potential model. 



length of the slowing region then would need to be about 
L w 3 (I) 2tN w 1 m, which is a reasonable length for a 
molecular beam apparatus. 

Alternatively, the additional lasers for repumping the 
X(v = 1,3,4 and/or 5) levels might be avoidable if 
other techniques for beam deceleration such as stimulated 
slowing [lU or Stark deceleration [l9| were employed to 
slow the beam. However, each of these would certainly 
introduce additional complications. Overall, it is diffi- 
cult to estimate the possible improvement of sensitivity 
to T-odd effects with an approach like this, since it will 
depend on details such as efficiency of trap loading, life- 
time in the trap, etc. Nominally, to achieve a factor of 
^100 improvement in sensitivity would require trapping 
~ 10 7 molecules with lifetime ~ 1 sec. This is similar to 
what has been achieved using Stark deceleration to load 
a storage ring of molecules, even without the additional 
advantage of laser cooling (albeit using species with pa- 



rameters chosen to optimize the efficiency of slowing and 
trapping) [201 ]. Hence, we believe that such an approach 
may be viable and deserves further investigation. 

In summary, the X 1 Yi + — £? 3 ni transition looks quite 
promising for the laser cooling of TIF. While the simulta- 
neous operation of several narrowband ultraviolet lasers 
will be costly, present laser technology can provide ade- 
quate power at the required excitation frequencies. If the 
cooling is successful, this should open up the possibility 
of a new generation of fundamental symmetries experi- 
ments in TIF with significantly improved sensitivity. 
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